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Abstract. We present selected results obtained by the ALICE Collaboration concerning 
charmed hadron and quarkonium production in Pb-Pb collisions at the LHC. 
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1. Introduction 

The goal of high-energy nucleus- nucleus collisions is to produce and characterize a state of nuclear 
(QCD) matter at (energy) densities well above the nuclear ground state (eq —0.15 GeV/fm 3 ). At 
high densities and/or at high temperatures one expects [TJ [2] that quarks are no longer confined 

.ys in protons and neutrons but move freely over distances larger than the size of the nucleon (~1 

fm=10~ 15 m). Employing Quantum Chromo-Dynamics calculations on a space-time lattice, a 

^& deconfinement phase transition for an energy density of about 1 GeV/fm 3 was predicted (see 

[3j for an early review). The characterization of quark-gluon matter in terms of its equation of 
state (EoS, relating pressure to energy) and of its transport properties (such as viscosity) and 
delineating its phase diagram is a major ongoing research effort in collisions of nuclei at high 
energies. 

Heavy quarks (charm and bottom) are a prominent part of the so-called "hard probes", 
observables characterized by an energy (mass) scale much larger than the temperature of 
the deconfined medium (of a few hundred MeV). These observables are produced in primary 
hard scattering processes and their yields are calculable in elementary hadronic collisions with 
perturbative QCD techniques [H [5j El |7]- The production of hard probes in nucleus-nucleus 

(A A) collisions is quantified by the nuclear modification factor: Raa = to — yd^jv P Jdv &v ' wnere 
d 2 N/dp<rdy denotes the transverse momentum (px) and rapidity (y) differential yield of a given 
observable measured in AA or pp collisions and (N co n) is the average number of nucleon-nucleon 
collisions over the given centrality interval of AA collisions; (iV co u) is calculated using the Glauber 
model |8|. Also calculated within this geometric model is the number of participating nucleons, 
(Afpart)) use d in the following for study of centrality dependence of -Raa- 

The pr spectra of hadrons carrying light quarks and their azimuthal distributions with respect 
to the reaction plane (defined by the beam axis and the impact parameter of the colliding nuclei) 
exhibit features of collective flow and are described well by hydrodynamical models (see a recent 
review in |9j). In non-central collisions, the initial approximately elliptic shape (in the transverse 
plane) of the overlap zone of the colliding nuclei leads, through the initial gradients of the energy 



density (or pressure), to anisotropic spatial (azimuthal angular) emission of hadrons. This is 
quantified by the second order (quadrupole) Fourier coefficient v<i = (cos(2</>)), where <j) is the 
azimuthal angle with respect to the reaction plane. 

The study of (deconfined) QCD matter has entered a new era in year 2010 with the advent of 
Pb-Pb collisions at the Large Hadron Collider (LHC), delivering the largest ever collision energy, 
y / SNN=2.76 TeV, more than a factor of 10 larger than previously available. We present selected 
results obtained by the ALICE Collaboration concerning charmed hadron and quarkonium 
production. The ALICE detector is described in [ID]. The results presented below correspond 
to minimum-bias Pb-Pb run of the year 2010 and to the high-statistics Pb-Pb run of the year 
2011. An integrated luminosity Li n t —30 fib^ 1 was collected in 2011 with centrality triggers for 
the ALICE Central Barrel at midrapidity (|y| < 0.9), while Li nt ~70 /xb" 1 was collected with 
(di)muon triggers in the forward Muon Spectrometer (2.5 < y < 4). 



2. Open charm hadrons 

Produced very early in the collision, at times t ~ l/2m q , heavy quarks experience the full 
evolution of the hot fireball of partons. At variance with light quarks, which can be thermally 
produced, heavy quarks maintain their identity through the hot stage of the collision. They help 
to address the question if parton energy loss (by gluon radiation) exhibits the expected quark 
mass pattern. The theoretical expectation is that heavy quarks (charm and bottom) lose less 
energy (by gluon radiation) compared to lighter (up, down, strange) quarks [IT]. The related 
question, often asked, is whether heavy quarks thermalize alongside the light quarks and gluons. 
The measurements at the LHC were eagerly awaited in order to establish in a new energy regime 
the observation at RHIC, at a center of mass energy per nucleon pair ^/snn = 200 GeV, of heavy 
quark energy loss and collective flow [12j . 
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Figure 1. Invariant mass distributions for charmed meson candidates for central (0-7.5% 
centrality) Pb-Pb collisions. Particles and antiparticles are included. 

The production of hadrons carrying open charm is studied in ALICE via the direct 
reconstruction of hadronic decays [13] at midrapidity, |y| < 0.5, and via the measurement 
of decay electrons [H] (at midrapidity) or muons [15] (at forward rapidity, 2.5 < y < 4). 
Invariant mass distributions for charmed meson candidates are shown in Fig. Q] for central (0- 
7.5% centrality) Pb-Pb collisions. Owing to the very good resolution of track impact parameter 
in the transverse plane (65 /im for p^ = 1 GeV/c and reaching asymptotically 20 y«m) and to the 
good hadron identification performance, a very good signal to background ratio S/B of about 
10% is reached for various D-meson species. The very good momentum resolution of ALICE at 
low px leads to invariant mass resolutions around 10 MeV/c 2 for the D + and D° mesons. 

The transverse momentum dependence of the nuclear modification factor Raa of D mesons 
in central Pb-Pb collisions is shown in Fig. [2] [14J. The left panel shows the measurements for 
three D-meson species, D + , D°, and D* + , while the right panel shows the average values of these 
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Figure 2. The transverse momentum dependence of the nuclear modification factor of prompt 
D mesons (particles and antiparticles) in central (0-7.5% centrality) Pb-Pb collisions (left panel). 
The right panel shows the average values of -Raa f° r D + > D°, and D* + meson species in 
comparison to data for all charged particles and for charged pions (for 0-10% centrality). 



species in comparison to data for all charged particles and for charged pions (which are for 0-10% 
centrality) . The similar Raa values of D mesons and of pions (or charged particles) , well below 
unity indicate that the energy loss of charm quarks is of similar magnitude compared to that 
of lighter quarks or gluons. The first measurement of Raa f° r D s mesons in heavy-ion collisions 
[16] is also shown in Fig. [2] (left panel). In the highest measured pp bin (8-12 GeV/c), the .Raa 
for D s mesons is compatible with that of non-strange charmed mesons. At lower pi, the .Raa 
of D s mesons seems to increase, but with the current statistical and systematic uncertainties no 
conclusion can be drawn on the expected enhancement of D s mesons with respect to non-strange 
D mesons |17| . 

The transverse momentum dependence of the elliptic flow coefficient V2 of charmed mesons 
|14] is shown in Fig. [3] together with that of inclusive charged hadrons for the centrality interval 
30-50%. Despite the large uncertainties, the challenging measurement of V2 for D mesons shows 
a significant anisotropy, suggesting that charm quarks may take part in the collective motion 
built up at the quark level in the deconfined stage. 

The -Raa and V2 results are confronted with model predictions in Fig. HJ The theoretical 
models [TSJ 113 1201 ED E21 1231 E3] implement parton energy loss via elastic (collisional) and 
radiative mechanisms. The models are able to describe the features of the data, but a 
quantitative description of both energy loss and elliptic flow remains challenging for all models. 
As noted in [23] , an increase of -Raa as a function of p-r is generic for a constant fractional parton 
energy loss and for a parton spectrum of power-law form, with the power exponent increasing 
as a function of px- 

A further challenge for theory is the consistent description of energy loss for light [25] and 
heavy quarks [13], illustrated in Fig. [5] (left panel). Partonic energy loss models achieve a 
good description at high p^ while the low p^ region is not described well. The data suggest 
-Raa > R\ai as expected. Further hints on the quark flavor dependence of the nuclear 
modification factor are given in Fig. [5] (right panel) for intermediate p^ values (where the .Raa 
values exhibit a shallow minimum, see Fig. 0]). The ALICE data on light flavors [25J and charmed 
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Figure 3. The transverse momentum dependence of the elliptic flow coefficient V2 of charmed 
mesons for the centrality interval 30-50%. The data for inclusive charged hadrons is included 
for comparison. 
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Figure 4. Transverse momentum dependence oi Raa (left panel, central collisions) and V2 (right 
panel, 30-50% centrality) of D mesons. The mesurements are compared to model predictions 
(see text). 



hadrons [H] are compared to the measurements of primordial 3/tp and J/ip from decays of B 
hadrons [26]. Although the details of such a comparison need to be carefully considered [23], the 
data suggest the expected hierarchy of flavor dependence of in-medium parton energy loss. A 
further challenge to theory is the extraction of the heavy quark (momentum) diffusion coefficient, 
a quantity which was recently calculated within lattice QCD [27] . 



ALI-PREL-336BS 



+1 < 

B < 

rr 
cc 



- Pb-Pb, \ls^,=2.76 TeV 
45 1. Centrality: 0-20% 

: Average D°,D*,D"*, |y|<0.5' 
4 



Rad (Vitev) 

Rad + dissac (Vitev) 

WHDG rad+coll (Horewitz) - 

AdS/CFT Drag (Horowitz) 

CUJETI.O(Buzzatti) 

I 




0.5 -_ 



8 10 



12 14 16 

P T (GeV/c) 




ALICE 
Pb-Pb, 4KZ = 2.76 TeV - 



• Average D, D + , D * 6<p <1 2 GeV/c 
Charged particles, 6<p <1 2 GeV/c 
^ HF muons, 2.5<y<4.0, 6<p <10 GeV/c 

CMS ncn-prompt JA|/, |y|<1 .2, p >6.5 GeV/c 
[|] CMS prompt J/y/, p >6.5 GeV/c 

(CMS read from CMS-PAS-HIN-1 2-01 4, norm. unc. 9% not shown) 



0, 



common normalization uncertainty on ALICE data: 7% (peripheral) to 4% (central]^ 

50 100 150 200 250 300 350 400 

<N ) 

part 



Figure 5. Left panel: transverse momentum dependence of the ratio of Raa of D mesons 
and pions, data and model predictions. Right panel: centrality dependence of -Raa for hadrons 
carrying light and heavy quarks (see text). 



3. Quarkonium 

Among the various suggested probes of deconfinement, charmonium plays a distinctive role. 
It is the first hadron for which a clear mechanism of suppression in the QGP was proposed 
early on, based on the color analogue of Debye screening |28| (see |29| for a recent review). A 
suppression compared to pp collisions was observed in pA collisions [30] , which was understood 
as a destruction of the pre-resonant cc state by the nucleons of the colliding nuclei. The 
measurements in Pb-Pb at the SPS (^/snn = 17-3 GeV) [31J and in Au-Au at RHIC j3"2] 
demonstrated an "anomalous" suppression, attributed to the dense and hot QCD matter. The 
theoretical interpretation was a sequential suppression [33 ,33] of various charmonium states as a 
function of energy density (temperature) of the QGP. In this picture, the J ftp meson survives in 
the QGP and only the contribution from the ip(2S) and \ c charmonium states (leading to about 
10% and 30% of the total J/tp yield, respectively) vanishes as a consequence of the melting of 
ip(2S) and Xd,2 states. As appealing as this interpretation is, it was questioned by the non-zero 
ip(2S) yield measured at the SPS [35J and also later by the re-evaluation of cold nuclear matter 
suppression derived from p-A measurements |30j . 

The ip(2S) yield relative to the yield of J/ip was noted earlier |36[ [37] to correspond, for 
central collisions, to the chemical freeze-out temperature derived from fits to other hadron 
abundances. This led to the idea of statistical hadronization of charm quarks in nucleus- nucleus 
collisions [37J. In this model, the charm quarks produced in initial hard collisions thermalize in 
the QGP and are "distributed" into hadrons at chemical freeze-out. All charmonium states are 
assumed to be fully melted (or, more precisely, not formed at all) in the QGP and produced, 
together with all other hadrons, exclusively at chemical freeze-out. This model [38(. I39| (see also 
a recent review [30]) has gained support from experimental data at RHIC [32^ 141] . The model 
predicted a notably different pattern of 3/ip production at the LHC. Depending on the charm 
production cross section, even an enhanced production relative to pp collisions could be expected 
at the LHC in central Pb-Pb collisions [39] (see [12] for further predictions). The statistical 
model predicted [13] an energy-independent relative production ratio ip{2S) / i / ip about 4 times 
smaller in (central) AA collisions compared to pp. 

Proposed at the same time as the statistical hadronization model, the idea of kinetic 
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recombination of charm and anti-charm quarks in the QGP |44J is an alternative quarkonium 
production mechanism. In this model (see [451 146] for recent results), a continuous dissociation 
and regeneration of charmonium takes place in the QGP over its entire lifetime. Besides 
the charm production cross section, this model has as input parameters the time dependence 
of the temperature of the fireball as well as relevant cross sections and assumptions on the 
melting scenarios of charmonium states. Important observables, like the transverse momentum 
dependence of production yields and of elliptic flow can be calculated within the transport 
(kinetic) models. It predicted a rather small J/tjj regeneration component at RHIC and a 
sizable one at the LHC. 

The measurement of J/ip production in Pb-Pb collisions at the LHC was expected to provide 
a definitive answer on the question of (re)generation. The first data, measured at high-pr by 
ATLAS [47] and CMS [4S], showed a pronounced suppression with respect to pp collisions. 
Subsequently, it was seen that this J/tp suppression is of the same magnitude as that of 
open-charm hadrons [13] , as shown above (Fig. [41 right panel). This may indicate that the 
high-pT charm quarks that form either D or J/?/; mesons had the same dynamics including a 
thermalization stage and a late hadronization. 
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Figure 6. Invariant mass distributions measured for the /i + fi~ channel (Muon Spectrometer, 
2.5< y <4, right panel; the dashed line is a fit to the background) and e + e~ channel (Central 
Barrel, \y\ <0.9, right panel; the lower panel is the signal obtained after the mixed-event 
background subtraction, with the line the signal form obtained from Monte Carlo simulations). 



Quarkonium is measured in ALICE via its dielectron channel at midrapidity and via the 
dimuon channel at forward rapidity, respectively, as shown in Fig. [6] (see [49 [ [50] for the details 
of the analysis). A clear difference to the RHIC results [32[ I4T] was seen in the first LHC 
measurement of the overall (inclusive in p^) production, performed in ALICE [49J. In this case, 
where low px J/V production is dominant, less nuclear suppression (larger Raa value) is seen at 
the LHC both at forward rapidity [51J and at mid-rapidity [52J, as shown in Fig.[7J Corroborated 
with the CMS data [48] , teh result indicates an enhanced production of J/tp at 1ow-j?t compared 
to high px- 
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At the LHC, the estimated energy density is at least a factor of 3 larger than at RHIC 
[53], leading to an initial temperature most likely above the one required for 3/tjj dissociation. 
Therefore, one can conclude that the production mechanism of 3/ip and charmonium in general 
at the LHC is determined (to a large extent) by regeneration in QGP or by generation at 
chemical freeze-out. 
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Figure 7. Centrality dependence of the nuclear modification factor for inclusive 3/ip production. 
The ALICE measurement (preliminary) at the LHC is compared to the PHENIX data at RHIC. 
The two panels show the data at forward rapidity (left) and at mid-rapidity (right). 
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Figure 8. Centrality dependence of the nuclear modification factor for 3/ip at all momenta 
in comparison to theoretical models, for forward rapidity (left panel) and mid-rapidity (right 
panel) . 



Indeed, both the statistical hadronization |42j and transport |45} I46] models reproduce the 
data [49 1, as seen in Fig. [HJ Based on these observations, the 3 ftp production can be considered 
a probe of QGP as initially proposed [28], but may not be a "thermometer" of the medium 
[34] . Within the statistical model, the charmonium states become probes of the phase boundary 
between QGP and hadron phase. This extends with a heavy quark the family of quarks employed 
for the determination of the hadronization temperature (via the conjectured connection to the 
chemical freeze-out temperature extracted from fits of statistical model calculations to hadron 
abundances). 
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Figure 9. . Transverse momentum dependence of the nuclear modification factor for J/ip at 
forward rapidity for the centrality range 0-20%. The preliminary ALICE data at the LHC are 
compared to measurements at RHIC by PHENIX. 



The transverse momentum dependence of the nuclear modification factor, shown in Fig. [9j 
is, at the LHC, dramatically different than the one measured at RHIC. At low-px the nuclear 
suppression is significantly reduced (i.e. larger Raa) at the LHC |51j compared to RHIC |32j. 
Transport model calculations reproduce the data quantitatively, as can be seen in Fig. [TO] (with 
model of ref. [46J; the open bands represent the yield due to regeneration). In current models 
[33 SSj, about half of the low-pp J/ip yield in Pb-Pb collisions at y / SNN = 2.76 TeV is produced 
by the recombination of charm quarks in QGP, while the rest is due to primordial production. 
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Figure 10. Nuclear modification factor (left panel, for centrality ranges 0-20% and 40-90%) and 
elliptic flow (right panel, for centrality range 20-60%) of J/ijj mesons as a function of transverse 
momentum. The ALICE data is compared with parton transport model predictions. 



Both the kinetic and the statistical hadronization models require thermalization of the charm 
quarks in QGP. As a consequence they will follow the collective behavior of the bulk QGP and 
their flow will be reflected in that of charmed hadrons and quarkonia. Indeed, elliptic flow of 
J/tfi at LHC energies was predicted within a transport model [15]. The first measurement at 
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the LHC [54], shown in Fig. [10] (right), provides a tantalizing hint of a non-zero elliptic flow of 
J/ip (the significance of non-zero flow for the centrality range 20-60% is 2.3a in the pt range 
2-4 GeV/c). The preliminary ALICE data is consistent with the expectation from transport 
models, but precision data is needed in order to be able to extract information on the QGP 
properties and on the amount of J/ip produced via regeneration. 

The picture outlined above for the charmonium production, extracted by comparison of data 
and model predictions, remains to be tested by precision measurements at the LHC, which will 
place further constraints to models. A generic feature expected within (re)generation models is 
an increase of Raa at the higher energy (^/snn = 5.1 TeV) of collisions expected in the LHC 
after the current shutdown period. In particular, both the statistical and the transport models 
employ as input parameter the cc production cross section <t c5 . The sensitivity of the calculations 
to a C c is rather large (see Fig. [5]). A precision measurement of a C c in Pb-Pb collisions, within 
reach with the proposed ALICE upgrade [55], will place an important constraint to models. A 
precision measurement of excited states of charmonium, in particular of ip(2S), for which current 
data |56] are rather imprecise, could allow to discriminate between the regeneration [46] and 
statistical production [43] mechanisms. 

4. Summary and outlook 

A wealth of data has been measured by ALICE on charmed hadrons and charmonium in Pb-Pb 
collisions at the LHC. Parton energy loss is strong also for charm quarks and its quantitative 
description together with the observed elliptic flow in theoretical models remains a challenge. 
Charmonium exhibits features of production via regeneration in the hot deconfined medium or 
via statistical hadronization at the phase boundary. Disentangling these two pictures is the next 
challenge for the experiments and theory. 
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